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EFFECTS  OF  MAGNETIC  AND  COLLISIONAL  VISCOSITY 
ON  THE  EXB  GRADIENT  DRIFT  INSTABILITY 

I.  INTRODUCTION 

Much  has  been  written  in  the  recent  past  regarding  the  'freezing' 
phenomenon  exhibited  by  ionospheric  ion  clouds.  Recently  the  debate  has 
centered  upon  whether  the  explanation  for  the  phenomenon  is  a  2-  or  3- 
dimensional  effect.  We  shall  not  enter  that  debate  here,  but  rather 
investigate  one  of  the  more  promising  2-D  effects,  proposed  by  Sperling 
and  Glassman  (1985).  They  propose  that  collisional  and  magnetic 
viscosities,  and  respectively,  [Braginskii,  1965]  suppress  the 

growth  of  the  gradient  drift  instability  at  short  wavelengths,  and  that 
this  effect  can  explain  the  'freezing'  phenomenon  in  both  barium  ion 

clouds  and  nuclear  plumes.  The  nonlocal  stability  analysis  based  on  the 
full  stress  tensor  presented  in  Sperling  and  Glassman  (1985)  considered 

only  the  asymptotic  case  of  kL  <<  1  where  k  is  a  Fourier  wavenumber  and  L 
is  the  gradient  scale  length.  The  results  presented  in  their  paper  were 
quite  intriguing  and  showed  that  and  can  exhibit  a  strong 

stabilizing  effect  on  the  gradient  drift  instability. 

However,  one  must  seriously  question  the  applicability  of  the  above 
work  to  barium  clouds  or  nuclear  striations  since  the  results  are  valid 
only  in  the  long  wavelength  limit  which  corresponds  to  L  ->  0,  and  the 

stress  tensor  itself  is  valid  only  for  L  >>  where  pi  is  the  ion 

gyroradius.  What  are  the  instability  characteristics  for  finite  L,  and  in 
particular  for  kL  -  0(1)  as  we  expect  it  to  be  for  some  cases  of  interest. 
This  poses  a  difficult  numerical  problem,  and  has  not  henceforth  yielded 
to  solution.  However,  we  have  recently  developed  numerical  techniques 
which  allow  us  to  solve  the  linearized  equations  for  the  full  stress 
tensor,  and  to  investigate  the  applicability  and  validity  of  the  long 
wavelength  asymptotic  results.  We  find  that  the  gradient  scale  length  L 
does  have  a  strong  effect  on  the  growth  rate,  and  thus  that  long 
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wavelength  analysis  is  often  inadequate  for  realistic  parameters. 


II .  FUNDAMENTAL  EQUATIONS 

The  equations  describing  a  two  dimensional  plasma  in  a  plane 
perpendicular  to  an  ambient  magnetic  field  B  have  been  given  in  many 
places  [e.g.,  Mitchell  et  al.,  1985].  What  is  new  here  is  the  presence  of 
terms  which  account  for  magnetic  and  collisional  viscosity.  As  with 
Sperling  and  Glassman  (1985),  we  shall  ignore  all  electron  collisions.  We 
shall  do  this  because  we  do  not  wish  to  address  the  effects  of  these 
collisions  (primarily  diffusion)  here,  and  because  the  diffusion  induced 
by  these  collisions  precludes  the  attainment  of  a  true  equilibrium  upon 
which  to  perform  a  stability  analysis. 

Consider  a  two  dimensional  plasma  consisting  of  ions  and  electrons, 
embedded  in  a  neutral  gas  and  in  a  constant  ambient  magnetic  field  B 
perpendicular  to  the  plasma  plane.  If  we  neglect  ion-electron  collisions, 
the  momentum  equation  describing  species  a  is 


v  x  B 


fir  +  v  •  y)v  -  2S  (E  ♦  =SLUi)  .  v  fv  -  ul  -  -IS-  ♦  5 
(3t  ~  a  )  ~a  m  c  )  an  (~a  ~n)  n  m  a 


+  if  <l> 

a 


where  the  subscript  a  denotes  the  plasma  species  (i  for  ions,  e  for 
electrons,  for  example),  n  is  the  species  number  density,  v  is  the  species 
fluid  velocity,  P  *  n  k_T  is  pressure,  E  is  the  electric  field,  g  is  the 
gravitational  acceleration,  q  is  the  species  charge,  is  the  species 
collision  frequency  with  the  neutral  gas,  Un  is  the  neutral  wind  velocity, 
c  is  the  speed  of  light,  kg  is  Boltzmann's  constant,  m  is  the  species 
particle  mass,  and  L  is  the  "force"  per  particle  due  to  the  and 
terms  in  the  Braginskii  stress  tensor.  We  can  rewrite  this  equation  as 
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If  we  place  ourselves  in  a  Cartesian  coordinate  system  in  which  B  is 
aligned  along  the  z  axis  then  (2)  yields 


F  x  e  . 
-a  z 


(6) 


Strictly  speaking,  this  equation  only  applies  to  the  perpendicular 
component  of  v  .  We  assume  F^  =0  so  that  v^  =  0.  Note  that  (6)  above 
is  actually  an  implicit  expression  for  v^,  since  F^  is  a  function  of  v  . 
Hence  v^  is  solved  iteratively.  We  define 


and 


v  ,  =  — =■  F  t  x  e 
-oil  q  B  -od  z 


~a2  =  ft  -o^^al5  X  Sz 


Then  can  be  approximated  as 


V  =  V  .  +  V  _ 

-a  ~al  ~<x2 


(7) 


(8) 


(9) 


The  success  of  this  procedure  obviously  depends  on  F^  being  small  with 
respect  to  F^,  i.e.,  v^  <<  v^. 
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The  reader  vill  perhaps  be  disturbed  by  the  absence  of  explicit 
collision-frequency-dependent  Pedersen  and  Hall  mobilities.  Indeed,  the 
above  procedure  is  not  the  one  usually  followed,  wherein  the  term  v  m  v 
is  included  in  (2)  rather  than  (5)  (e.g.,  Zalesak  et  al.,  1985).  The 

procedure  used  here  is  less  accurate  than  this  usual  approach  (e.g.,  Hall 

currents  due  to  electric  fields  are  totally  absent  here).  However,  we  use 
this  simplified  set  here  for  three  reasons.  First,  it  is  accurate  as  long 
as  v^/C^  <<  1  and  |^2  I  ^  second,  it  is  simple;  and  third,  it  is 

the  approach  used  by  Sperling  and  Glassman  (1985),  with  whom  we  wish  to 
compare  our  results. 

We  now  list  the  additional  assumptions  we  need  to  make  in  order  to 

recover  the  linearized  equations  of  Sperling  and  Glassman  (1985).  Ue 

assume  that  both  the  ion  temperature  T^  and  the  electron  temperature 
are  constants  in  space  and  time,  neglect  all  electron  collisions,  neglect 
gravity  g,  assume  quasi-neutrality  (  ni  =  n  *  n  ),  assume  singly  changed 
ions  (  q^  *  -  qg  =  e  ),  and  assume  that  the  electric  field  is 
electrostatic  (i.e.,  E  =  -  7$) .  The  smallness  of  mg,  together  with  our 
neglect  of  electron  collisions,  allows  us  to  neglect  all  but  the  first  and 
third  terms  on  the  right-hand  side  of  (1)  for  electrons.  Thus 


v 

~e 


=  v^i  =  |  [E  +  (kBTe/e)(Vn/n)]  x  ez  (10) 


since  F  g  =  0. 


The  ion  velocities  are  given  by 


~i2  =  eB 
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To  solve  for  the  evolution  of  the  plasma,  ve  shall  need  tvo  equations: 
(1)  a  continuity  equation  for  electrons  or  ions,  and  (2)  an  equation  for 
current  continuity  and  quasi-neutrality,  i.e.,  V  •  J  =  0  where  J  is  the 
electric  current  en(v^  -  v^).  We  choose  the  electron  continuity  equation 
for  simplicity: 


3n 

3t 


(nv^)  = 


( c  ‘ 

7  •  r-  n7$  x  e 

z) 


(18) 


The  current  density  is  given  by 
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(19) 


thus 
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7  hti2  *  *  *z)]  <20> 

If  we  denote  the  x  and  y  components  of  as  a  and  v,  respectively,  and 
denote  partial  differentiation  by  subscripts,  then  we  can  write 


L(v..)  =  - 

~v  il'  n 


(r  D)x  *  Hi  *  Mx  *  My  -  M; 
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where 
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Here  is  the  ion-ion  collision  frequency,  A  is  the  ion  atomic  mass, 
and  is  the  ion  cyclotron  frequency  [Braginskii,  1965,  Sperling  and 
Glassman,  1985 ] .  T\  in  (31)  is  expressed  in  eV. 

Noting  that 


*il  *  * 


„  v.  m.  . 

e  =  |  (7*  -  — — -  U  ) 
z  B  l  e  -nj 


(32) 


we  obtain  equations 
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-  mi(  It 
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where 
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(39) 


The  final  equations  describing  our  plasma  are  then 


It  ~  I  (7*x  •«)*  7n  =  0  (40) 


]v* 
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~Q 
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2  2 

Note  that  we  have  dropped  the  small  term  V  •  [ en Un n / ]  in  (41)  to  be 
consistent  with  Sperling  and  Glassman  (1985).  Note  also  that  in  the  above 
equations,  electron  and  ion  temperature  can  change  the  plasma  evolution 
only  through  and  q^.  This  is  a  result  of  our  assumptions  of  no 
electron  collisions,  and  of  spatially  uniform  ion  and  electron 
temperatures. 
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III. 


STABILITY  ANALYSIS 


We  consider  the  following  equilibrium  condition,  which  we  denote  by 
zero  subscripts.  We  take  ng,  and  hence  ri^Q  and  to  be  functions  of  y 
only,  and  Un  to  be  constant  and  in  the  y  direction 


5n0  5n10  _  3n30 

3x  ~  3x  ~  3x 


(42) 


U 

-n 


U 

n 


e 


y 


This  equilibrium  configuration  is  characterized  by 


* 0  =  0 


(43) 


(44) 


v. 

m  ,, 

v  -  in  =  n -  ^  e 

-lio  2.  n  x 


(45) 


(The  reader  may  have  noticed  that  this  appears  not  to  be  an  equilibrium 
since  v^q  *  this  is  a  consequence  of  errors  introduced  in  the  Sperling 
and  Glassman  (1985)  ordering.  If  the  treatment  had  been  the  more  accurate 
classical  one,  a  true  equilibrium  would  be  obtained.) 

We  introduce  the  perturbation  quantities 


Linearizing  (40) 


n(x,y, t)  =  nQ(y)  +  n(y' t)e 

<Kx,y,t)  =  <ky,t)e 

and  (41)  we  find  that 
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The  usual  approach  to  stability  analysis  is  to  use  the  following 
substitutions 


| — »  ik. 
3x 


(53) 


37  *  1W 


(54) 


in  (48)  and  (49);  thereby  eliminating  either  n  or  and  obtaining  an 
eigenvalue  problem  for  a  fourth  order  differential  operator.  This  is  the 
approach  used  in  Sperling  and  Glassman  (1985).  which  yielded  an  equation 
they  solved  only  in  the  long  wavelength  limit.  Here  we  solve  the  problem 
with  no  approximations  by  treating  equations  (48)  -  (52)  as  an  initial 
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value  problem,  using  a  random  seed  perturbation.  All  derivatives, 
including  temporal  derivatives,  are  discretized  using  standard  finite 
differences.  Assuming  that  there  is  a  fastest-growing  eigenmode,  this 
mode  will  emerge  from  the  noise,  eventually  reaching  a  point  where,  on  a 
relative  scale,  it  is  virtually  the  only  mode  left  in  the  system.  At  this 
point  both  the  fastest  growing  eigenmode  and  its  complex  eigenf requency 
have  been  isolated.  Before  we  move  on  to  these  results,  let  us  look  at 
two  limiting  cases  for  which  solutions  are  somewhat  easier  to  obtain:  the 
asymptotic  cases  of  long-  and  short-wavelength  limits. 

IV.  ASYMPTOTIC  RESULTS:  SHORT-  AND  LONG-VAVELENGTH  LIMITS 
A.  Short-Wavelength  Limit 

We  consider  first  the  short-wavelength  limit  kL  >>  1,  where  L  a 

ng(3nQ/3y)~*  is  the  gradient  scale  length.  We  take  the  perturbation 
quantities  in  (48)  -  (52)  to  be  invariant  in  y.  Making  the  substitutions 
(53)  and  (54),  (48)  and  (49)  become 

3nn 

iai"  =  ~  |  ik  3^T  *  (55) 


pQ 


.2  -  T  qp0  ,2f.  Vin  ..  )  r  Vin  „  ..  - 

k  *  -  I  k  +  ik  5—  UJ*  =  —  Uneikn 


in 


nlOyyk2  +  n10k 


4  +  2r|30yik3 


(56) 


where  I  is  a  flag  on  ion  inertia  (I  =  1  to  retain  ion  inertia,  I  =  0  to 

neglect  it).  Solving  (56)  for  n,  and  substituting  into  (55)  we  obtain 
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'“[^iXek] 


-  "  1  ^  ^  (iw  +  ik  of  Un) 


'  r2  (n10yyk2  +  n10k4  +  2n30yik;3) 


_  £.  ik  !2o 

B  3y 


Noting  that  (2^/v^e)  a  „  =  (c/B)nQ,  dividing  through  by 


pO 

and  multiplying  by  iU  k,  we  obtain 

n 


id) 


-  fe2  -  (*“  *  Koyy 


*  1lOk  *  2n30ylk3] 


2  1  3n0 
nQ  3y  Un 


This  is  simply  a  quadratic  equation  for  w  of  the  form 

2 

Am  +  Bw  +  C  =  0 


with 


A  -  Iv 


-1 

in 


B  =  -  i  +  Iv. “1  k  giS  (J 
in  2^  n 


ng  3y  n 


-infr  -si  -if, 

2,.  Un  ~  1  02  ffp0  (n 


b2  up0  V'lOyy  +  niOk  +  2 


(57) 

this  quantity, 

k2 

(58) 

(59) 

(60) 

*bOy*k) 

(62) 
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An  instructive  look,  at  how  the  and  h2  terms  might  contribute  to 


the  stabilization  of  the  gradient  drift  instability  is  obtained  by  first 
dropping  ion  inertia  (setting  I  =  0),  and  looking  at  the  collisional 
limit.  Then 


i _!2o 

“  "  lUn  n0  3y 


ci  -1/ 

%  2  V  ln 


pO  V'lOyy  +  n10k 


2ri30ylk) 


-1 


(63) 


The  growth  rate  y  is  given  by  minus  the  imaginary  part  of  u>: 

1 
B" 


1_  ^0 
Un  nn  9y 


o2  ffp0  (n10yy  +  n10k  ) 


1  +  g2  V>(n10yy  +  ^lO1*1  ) 


-  c  -L  . 
2  _2  ffp0  n30yk 

D 


(64) 


One  can  note  that  the  h2  terms  always  act  to  stabilize  the  collisional 
short-wavelength  limit,  but  that  the  terms  need  not  always  be  similarly 
stabilizing. 

In  the  limit  of  small  v^/52.^  (low  electron  density)  (25)  -  (31)  yield 


„  ,  ,  nikBTi  Vii  0  3  nk  T  M 

h,  =  1.2  =  0.3  nkBTi  ^2 


(65) 


n.k-T.  v. .  . 

*3  -  T71  25 1  =  2  nkBVsi 

11  1 


Since  the  mean  thermal  ion  gyroradius  is  given  by 


k  T 
2  B  l 

°i  =  - 2 

1  m.  52. 

l  l 


and  since 


2  .  52.  n  2  52. 

c  -1  l  B  c  _  l  c _ 

_2  ffp  ”  v.  nee  n2  ”  \>.  Bne 
B  y  in  B  in 


(66) 


(67) 


(68) 
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then 


,  c2  -1  „  -  2i  c _  3  (1  2  m  _  'j  2i  2  1  3n0 

b2  CTpO  n30y  *  \>.n  BnQe  3y  (.2  Opi  i Qi J  =  v.n  pi  nQ  3y 


(69) 


Thus,  if  ve  let  v^/2^  •»  0  ve  get 


Y 


VL 


.  2.2  4 ..  2  2 

2i  k  pi  /L  vin 


(70) 


where  we  have  defined 


(71) 


Recall  that  this  expression  is  valid  for  the  collisional,  short 
wavelength  limit  in  the  low  electron  density  regime  *  0).  Huba 
and  Zalesak  (1985)  have  shown  that  the  long  wavelength  limit  (kL  <<  1)  can 
often  be  obtained  from  the  short  wavelength  limit  by  making  the  simple 
substitution  1/L  -»  k(M  -  1)/(M  +  1),  yielding 


r  = 

■  -  3 ;  1  “of1  *  (1 

1  -  1)2  .  4 ' 4I — 1 

1 .  J  k  »ij 

(72) 

P.  *  1 

1  \>. 
m 

pi 

(73) 

Here  M  is  the 

ratio 

of  the  electron 

densities  on 

either  side  of  the 

discontinui ty . 

This 

collisional  long 

wavelength  limit 

has  been  obtained 

rigorously  by  Sperling  and  Glassman  (1985).  It  is  obvious  that  effects 
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will  strongly  damp  the  instability  whenever 


k.Pi  > 


.  /n. , 

in  1 ) 


\  1/2 


(74) 


in  this  short  wavelength  collisional  limit. 

B.  Long-Uavelength  Limit 

By  the  long-wavelength  limit  (kL  <<  1)  we  mean  the  solution  to  the 

linearized  equations  (47)  -  (52)  with 

n<  ?  y  <  y0 

nQ(y)  =  (75) 

n>  *  Mn<  *  y  >  y0 

The  concept  of  the  long  wavelength  limit  within  the  context  of  the 
stability  of  ionospheric  ion  clouds  was  introduced  by  Huba  and  Zalesak 
(1983),  and  has  become  a  standard  analysis  tool  in  the  field.  The  logic 
behind  the  use  of  the  long  wavelength  limit  rather  than  the  older  short 
wavelength  limit  is  that  the  clouds  in  question  tend  to  steepen 
considerably  before  structuring,  making  the  analysis  of  a  very  steep 
density  jump  a  more  plausible  approximation  than  that  of  a  smooth 
exponential  profile.  One  of  the  primary  purposes  of  this  paper  is  to 
compare  the  long  wavelength  results  with  those  of  continuous,  more 
realistic  profiles,  in  order  to  ascertain  the  applicability  of  the  long 
wavelength  approximation  for  realistic  clouds. 

The  equations  describing  the  long-wavelength  limit  are  derived  in 
Sperling  and  Glassman  (1985),  and  we  shall  not  repeat  that  analysis  here. 
Suffice  it  to  say  that  the  equations  are  easier  to  solve  than  the  full  set 
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(48)  -  (52),  and  that  we  have  developed  software  which  reproduces  the 
Sperling  and  Glassman  (1985)  results,  for  comparison  purposes. 


V.  NUMERICAL  RESULTS 

A.  Long-Wavelength  Comparisons 

All  of  the  results  we  report  here  are  for  an  atomic  oxygen  plasma  (m. 
*  16.0  nip)»  -  0.1  eV,  and  use  the  following  profile  for  n^ 


n0 


(76) 


5  -3 

n^  «  10  cm  ,  M  =  1000,  and  y q  is  the  center  of  our  computational  domain 

in  all  of  our  problems.  Three  different  values  of  \>.  a~e  used,  meant  to 

span  a  vide  range  of  altitudes:  vin  =  10  sec-1  (160  km  altitude),  v^n  =  1 

sec-*  (250  km  altitude)  and  \>in  *  0.1  sec-1  (500  km  altitude).  The 

neutral  wind  velocity  U  was  taken  to  be  U  =  -  100  m/s  in  all  cases. 

n  n 

Since  our  primary  goal  in  this  paper  is  to  examine  the  applicability 
of  long  wavelength  asymptotic  theory  to  realistic  situations  with  finite 
gradients,  we  shall  compare  calculations  in  which  only  Lq  changes  in  (76). 
Note  that  the  long  wavelength  limit  corresponds  to  Lq  ->  0. 

Figure  1  shows  results  for  the  case  =0.1  sec-*.  Displayed  is  a 

plot  of  r  vs.  k  for  Lq  =  0,  20  m,  100  m,  250  m,  500  m  and  1000  m.  Figures 
2  and  3  show  results  for  the  cases  \^n  =  1  sec-*  and  10  sec-*, 

respectively.  The  solid  curve  is  for  Lq  =  0  m  and  is  based  on  the  exact 

asymptotic  dispersion  equation;  the  curves  labelled  A,  B,  C,  D,  and  E  are 
for  Lq  =  20  m,  100  m,  250  m,  500  m,  and  1000  m,  respectively,  and  are 
based  on  numerical  solutions  to  (48)  -  (52).  Missing  data  corresponds  to 
cases  where  convergence  to  an  unambiguous  fastest-growing  eigenmode  could 
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not  be  obtained  with  the  present  code. 

It  is  obvious  from  Fig.  1-3  that  the  value  of  L  does  substantially 

affect  the  growth  rate  y.  There  are  two  equally  important  aspects  of  this 

issue  which  bear  on  our  discussion  here.  First,  there  is  the  absolute 

effect  of  L  on  y,  all  other  parameters  being  held  fixed.  From  a  look  at 

Figs.  1-3,  the  reader  can  see  that  the  growth  rate  for  a  finite  L  is 

always  smaller  than  that  for  L  =  0,  and  that,  in  general,  increasing  L 

decreases  y,  all  other  parameters  being  held  fixed.  The  degree  to  which 

this  reduction  takes  place  does  seem  to  be  a  function  of  the  ion-neutral 

collision  frequency,  with  the  effect  being  largest  for  large  \^  ,  and 

smallest  for  small  \>.  . 

in 

The  second,  perhaps  more  important,  aspect  of  the  effects  of  finite 
L,  is  the  degree  to  which  L  affects  the  k  for  which  r(k)  maximizes,  k 

max 

Ve  shall  also  be  interested  in  the  degree  to  which  the  curve  y(k)  is 

sharply  peaked  about  kmax.  Both  of  these  issues  are  highly  relevant  to 

the  freezing  models  which  have  been  proposed  [Glassman  and  Sperling,  1987, 

Zalesak  et  al.,  1988].  Again  referring  to  Figs.  1-3,  the  reader  will  note 

that  the  effect  of  increasing  L  is  to  make  k  smaller.  Since  k  is  in 

fact  the  freezing  scale  for  many  freezing  models,  this  is  a  very  important 

effect.  The  effect  of  increasing  L  on  the  degree  to  which  the  curve  r(k) 

is  sharply  peaked  is  more  difficult  to  describe.  As  L  is  increased  from  0 

to  20  m.  the  curve  becomes  more  sharply  peaked,  which  is  a  desirable 

attribute  for  the  freezing  model  of  Zalesak  et  al.  (1988).  However,  it  is 

clear  that  as  L  is  further  increased,  the  curve  eventually  becomes  less 

sharply  peaked  with  increasing  L.  Thus  larger  values  of  L  may  make  some 

freezing  models  which  depend  on  r(k)  being  sharply  peaked  about  k  less 

max 


desirable. 


B.  Short-Wavelength  Comparisons 

Having  compared  true  growth  rates  with  long-wavelength  asymptotic 
results,  and  found  the  asymptotic  results  lacking  to  some  degree, 
especially  at  short  wavelengths,  it  is  natural  to  assume  that  perhaps  the 
short-wavelength  limit,  (58)  -  (62),  might  be  a  fairly  good  approximation 
in  this  regime.  However,  it  is  not  clear  how  to  apply  (58)  -  (62)  to  our 
profile  (76),  since  the  spatial  derivatives  of  n^,  ,  and  ru  are 
functions  of  y.  We  have  chosen  to  evaluate  y(y)  as  given  by  (58)  -  (62) 
at  every  point  in  y  of  the  profile  (76),  and  to  take 


Yswl  =  max  y(y) ,  -  ®  <  y  <  +  « 


(77) 


In  Fig.  4  we  compare  r^-y^  with  the  numerically  computed  solution  for 


-1, 


the  case  L  =  1000  m,  for  all  three  values  of  \>.  (0.1,  1.0,  10.0  s-  ). 
Note  that  there  is  excellent  agreement  for  kL  >  10,  but  that  there  can  be 
substantial  disagreement  for  kL  -  1.  In  Figs.  5-7,  we  show  the  same 
comparison  for  the  cases  L  =  500  m,  L  =  250  m,  and  L  =  100  m, 
respectively.  Again,  note  that  kL  must  be  considerably  larger  than  1  to 
obtain  good  agreement  between  the  curves,  and  that  there  are  substantial 
errors  in  Yg^y  for  kL  -  1. 


in 


0.1)  of  Fig.  4-7  that  YgWL  seems  to  be 


We  introduce  some  bit  of  mystery  here.  The  reader  may  note  by 
comparing  curve  D  (v 
asymptoting  to  a  y  independent  of  L  for  large  k  in  this  highly 
collisionless  case.  We  do  not  as  yet  have  an  explanation  for  this 

phenomenon. 
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C.  Relative  Effectiveness  of  and 

Finally  in  this  section  we  wish  to  give  the  reader  a  feeling  for  the 
relative  effectiveness  of  the  ri^  and  terms  on  the  growth  rate  of  y.  In 
Fig.  8  we  show  a  plot  of  y  vs.  k  for  the  case  Lq  =  20  m,  and  all  three 
values  of  (0.1,  1.0,  and  10.0  s-^),  but  with  ^  =  0.  In  Fig.  9 
we  show  an  identical  plot,  but  for  ^  =  0  (with  taking  on  its  actual 
values).  In  Fig.  10  we  show  a  plot  with  both  and  set  to  actual 
values.  Note  that  the  achievement  of  peaked  values  of  y(k)  came  only  when 
tl2  terms  were  taken  into  account,  making  it  the  primary  piece  of  relevant 
physics  for  freezing  models  based  on  a  kmax-  This  is  verified  by  Fig.  11, 
where  we  show  a  similar  plot  with  =0,  and  only  tu  terms  "turned  on". 
We  do  not  as  yet  have  an  explanation  for  the  turning  up  of  y(k)  for  = 
0.1  at  very  large  k  in  Fig.  11.  (We  find  that  making  L  smaller  makes  the 
turning  up  vanish,  in  agreement  with  the  true  L  =  0  result.) 

VI.  CONCLUSIONS 

We  have  examined  the  role  of  magnetic  and  collisional  viscosity  (the 
0^  and  terms,  respectively,  in  Braginskii's  stress  tensor)  on  the  E  x  B 
gradient  drift  instability,  paying  special  attention  to  the  degree  to 
which  the  gradient  scale  length  L  affects  curves  of  y(k).  In  particular, 
we  wished  to  know  the  size  of  the  error  introduced  by  assuming  L  to  be 
zero  (the  long-wavelength  limit).  We  conclude  that  in  addition  to 
affecting  y  in  an  absolute  sense,  increasing  L  to  a  finite  value  decreases 
kmax»  the  k  for  which  y(k)  maximizes,  and  also  affects  width  of  the  peak 
of  y(k)  at  kmax.  Both  of  these  results  have  large  effects  on  freezing 
models.  We  have  also  examined  the  applicability  of  the  old  short- 
wavelength  limit,  or  local  analysis,  and  found  that  significant  errors 
would  be  made  unless  kL  is  very  large  (-  10). 
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Fig.  1  Plot  of  y  vs.  k  for  the  case  of  n^  =  10^,  M  =  1000,  Un  =  -  100 
m/s,  using  a  hyperbolic  tangent  density  p  ofile  and  =  0.1  s 

Curves  A,  B,  C,  D,  and  E  refer  to  L  =  1000  m,  500  m,  250  m,  100  m, 
and  20  m,  respectively.  Missing  data  corresponds  to  cases  for 
which  our  code  could  not  converge  to  an  unambiguous  fastest- 
growing  eigenmode.  The  solid  unlabeled  curve  is  the  long- 
wavelength  limit  (L  =  0). 
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K  ( KM— l ) 


Fig.  4  Plots  of  r  vs.  k.  for  the  case  of  =  10^,  M  =  1000,  un  =  -  10C 

m/s,  using  a  hyperbolic  tangent  density  profile,  and  L  =  1000  m. 
Curves  A,  B,  and  C  refer  to  \^n  =  0.1,  1.0,  and  10.0  s  , 
respectively.  Curves  D,  E,  and  F  refer  again  to  \k ^  =  0.1,  1.0, 
and  10.0  s-*,  respectively,  but  are  generated  using  the  short- 
wavelength  asymptotic  limit  (77),  ( 58 )  —  ( 62 ) . 


K  (KM- l ) 


Fig.  8  Plots  of  r  vs.  k.  for  the  case  of  n<  =  10^,  M  =  1000,  =  -  IOC 

m/s,  using  a  hyperbolic  tangent  density  profile,  and  L  =  20  m. 
Both  the  collisional  viscosity  and  the  magnetic  viscosity  hi 

have  been  neglected.  Curves  A,  B,  and  C  refer  to  \>^n  =  0.1,  1.0, 
and  10.0  s-^,  respectively. 
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OFFICER 

UNDER  SECY  OF  DEFENSE  FOR 
RESEARCH  AND  ENGINEERING 
DEPARTMENT  OF  DEFENSE 
WASHINGTON,  DC  20301 

01CY  ATTN  STRATETIC  &  SPACE 

SYSTEMS  (OS) 

COMMANDER/DIRECTOR 
ATMOSPHERIC  SCIENCES  LABORATORY 
U.S.  ARMY  ELECTRONICS  COMMAND 
WHITE  SANDS  MISSILE  RANGE,  NM  88002 
01CY  ATTN  DELAS-EO/F.  NILES 

DIRECTOR 

BMD  ADVANCED  TECH  CENTER 
HUNTSVILLE  OFFICE 
P.O.  BOX  1500 
HUNTSVILLE,  AL  35807 

01CY  ATTN  ATC-T/MELVIN  CAPPS 
01CY  ATTN  ATC-O/W.  DAVIES 

01CY  ATTN  ATC-R/DON  RUSS 


33 


PROGRAM  MANAGER 
BMD  PROGRAM  OFFICE 
5001  EISENHOWER  AVENUE 
ALEXANDRIA,  VA  22333 

01CY  ATTN  DACS-BMT/J .  SHEA 

COMMANDER 

U.S.  ARMY  COMM-ELEC  ENGINEERING 
INSTALLATION  AGENCY 
FT.  HUACHUCA,  AZ  85613 

01CY  ATTN  CCC-EMEO/GEORGE  LANE 

COMMANDER 

U.S.  ARMY  FOREIGN  SCIENCE  &  TECH  CTR 
220  7TH  STREET,  N.E. 

CHARLOTTESVILLE,  VA  22901 
01CY  ATTN  DRXST-SD 

COMMANDER 

U.S.  ARMY  MATERIAL  DEV  &  READINESS 
COMMAND 

5001  EISENHOWER  AVENUE 
ALEXANDRIA,  VA  22333 

01CY  ATTN  DRCLDC/J. A.  BENDER 

COMMANDER 

U.S.  ARMY  NUCLEAR  AND  CHEMICAL  AGENCY 
7500  BACKLICK  ROAD 
BLDG  2073 

SPRINGFIELD,  VA  22150 
01CY  ATTN  LIBRARY 

DIRECTOR 

U.S.  ARMY  BALLISTIC  RESEARCH 
LABORATORY 

ABERDEEN  PROVING  GROUND,  MD  21005 
01CY  ATTN  TECH  LIBRARY/ 

EDUARD  BAICY 

COMMANDER 

U.S.  ARMY  SATCOM  AGENCY 
FT.  MONMOUTH,  NJ  07703 

01CY  ATTN  DOCUMENT  CONTROL 

COMMANDER 

U.S.  ARMY  MISSILE  INTELLIGENCE  AGENCY 
REDSTONE  ARSENAL,  AL  35809 
01CY  ATTN  JIM  GAMBLE 

DIRECTOR 

U.S.  ARMY  TRADOC  SYSTEMS  ANALYSIS 
ACTIVITY 

WHITE  SANDS  MISSILE  RANGE,  NM  88002 
01CY  ATTN  ATAA-SA 

01CY  ATTN  TCC/F.  PAYAN,  JR. 

01CY  ATTN  ATTA-TAC/LTC  J.  HESSE 


COMMANDER 

NAVAL  ELECTRONIC  SYSTEMS  COMMAND 


WASHINGTON,  DC 

20360 

01CY 

ATTN 

NAVALEX  034/T.  HUGHES 

01CY 

ATTN 

PME  117 

01CY 

ATTN 

PME  117-T 

01CY 

ATTN 

CODE  5011 

COMMANDING  OFFICER 
NAVAL  INTELLIGENCE  SUPPORT  CENTER 
4301  SUITLAND  ROAD,  BLDG.  5 
WASHINGTON,  DC  20390 

01CY  ATTN  MR.  DUBBIN/STIC  12 

01CY  ATTN  NISC-50 

01CY  ATTN  CODE  5404/ J.  GALET 

COMMANDER 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO,  CA  92152 

01CY  ATTN  J.  FERGUSON 


NAVAL  RESEARCH  LABORATORY 


WASHINGTON,  DC 

20375- 

-5000 

26CY 

ATTN 

CODE 

4700/S. 

OSSAKOU 

50CY 

ATTN 

CODE 

4780/J . 

HUBA 

01CY 

ATTN 

CODE 

4701 

01CY 

ATTN 

CODE 

7500 

01CY 

ATTN 

CODE 

7550 

01CY 

ATTN 

CODE 

7580 

01CY 

ATTN 

CODE 

7551 

01CY 

ATTN 

CODE 

7555 

01CY 

ATTN 

CODE 

4730/E. 

MCLEAN 

01CY 

ATTN 

CODE 

4752 

01CY 

ATTN 

CODE 

4730/B. 

RIPIN 

2*CY 

ATTN 

CODE 

2628 

01CY 

ATTN 

CODE 

1004/P. 

MANGE 

01CY 

COMMANDER 

ATTN 

CODE 

8344/M. 

KAPLAN 

NAVAL  SPACE  SURVEILLANCE  SYSTEM 
DAHLGREN,  VA  22448 

01CY  ATTN  CAPT.  J.H.  BURTON 

OFFICER-IN-CHARGE 
NAVAL  SURFACE  WEAPONS  CENTER 
WHITE  OAK,  SILVER  SPRING,  MD  20910 
01CY  ATTN  CODE  F31 

DIRECTOR 

STRATEGIC  SYSTEMS  PROJECT  OFFICE 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  DC  20376 
01CY  ATTN  NSP-2141 

01CY  ATTN  NSSP-2722/ 

FRED  WIMBERLY 
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OFFICER  OF  NAVAL  RESEARCH 
ARLINGTON,  V A  22217 

OICY  ATTN  CODE  465 

OICY  ATTN  CODE  461 

OICY  ATTN  CODE  402 

OICY  ATTN  CODE  420 

OICY  ATTN  CODE  421 

COMMANDER 

AEROSPACE  DEFENSE  COMMAND/ XPD 
DEPARTMENT  OF  THE  AIR  FORCE 
ENT  AFB.  CO  80912 


OICY 

ATTN 

XPDQQ 

OICY 

ATTN 

XP 

AIR  FORCE 

GEOPHYSICS  LABORATORY 

HANSCOM  AFB,  MA 

01731 

OICY 

ATTN 

OPR/HAROLD  GARDNER 

OICY 

ATTN 

LKB/ 

KENNETH  S.W.  CHAMPION 

OICY 

ATTN 

OPR/ ALVA  T.  STAIR 

OICY 

ATTN 

PHD/ JURGEN  BUCHAU 

OICY 

ATTN 

PHD/ JOHN  P.  MULLEN 

AF  WEAPONS  LABORATORY 
KIRTLAND  AFB,  NM  87117 
OICY  ATTN  SUL 

OICY  ATTN  CA/ ARTHUR  H.  GUENTHER 

AFTAC 

PATRICK  AFB,  FL  32925 
OICY  ATTN  TN 

WRIGHT  AERONAUTICAL  LABORATORIES 
VRIGHT-PATTERSON  AFB,  OH  45433-6543 
OICY  ATTN  AAAI/VADE  HUNT 

OICY  ATTN  AAAI /ALLEN  JOHNSON 

DEPUTY  CHIEF  OF  STAFF 
RESEARCH,  DEVELOPMENT,  AND  ACQ 
DEPARTMENT  OF  THE  AIR  FORCE 
WASHINGTON,  DC  20330 
OICY  ATTN  AFRDQ 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01731-5000 
OICY  ATTN  J.  DEAS 

ESD/SCD-4 

COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WRIGHT- PATTERSON  AFB,  OH  45433 
OICY  ATTN  NICD/ LIBRARY 

OICY  ATTN  ETDP/B.  BALLARD 


COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
GRIFFIN  AFB,  NY  13441 

OICY  ATTN  DOC  LIBRARY/TSLD 

OICY  ATTN  OCSE/V .  CO TOE 

STRATEGIC  AIR  COMMAND/XPFS 
OFFUTT  AFB,  NB  68113 
OICY  ATTN  XPFS 

SAMSO/MN 

NORTON  AFB,  CA  02409 
(MINUTEMAN) 

OICY  ATTN  MNNL 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  MA  01731 

OICY  ATTN  EEP/A.  LORENTZEN 

DEPARTMENT  OF  ENERGY 
LIBRARY,  ROOM  G-042 
WASHINGTON,  DC  20545 

OICY  ATTN  DOC  CON  FOR 

A.  LABOUITZ 

DEPARTMENT  OF  ENERGY 
ALBUQUERQUE  OPERATIONS  OFFICE 
P.O.  BOX  5400 
ALBUQUERQUE,  NM  87115 

OICY  ATTN  DOC  CON  FOR 

D.  SHERWOOD 


EG AG,  INC. 

LOS  ALAMOS  DIVISION 

P.O.  BOX  809 

LOS  ALAMOS,  NM  85544 

OICY  ATTN  DOC  CON  FOR 
J.  BREEDLOVE 
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UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 
LIVERMORE,  CA  94550 


OlCY 

ATTN 

DOC  CON  FOR 
TECH  INFO  DEPT 

01CY 

ATTN 

DOC  CON  FOR 
L-389/R.  OTT 

01CY 

ATTN 

DOC  CON  FOR 
L-31/R.  HAGER 

LOS  ALAMOS  NATIONAL  LABORATORY 
P.O.  BOX  1663 

LOS  ALAMOS,  NM 

87545 

01CY 

ATTN 

J.  WOLCOTT 

OlCY 

ATTN 

R.F.  TASCHEK 

01CY 

ATTN 

E.  JONES 

01CY 

ATTN 

J.  MALIK 

01CY 

ATTN 

R.  JEFFRIES 

01CY 

ATTN 

J.  ZINN 

01CY 

ATTN 

D.  WESTER VELT 

01CY 

ATTN 

D.  SAPPENFIELD 

NATL.  OCEANIC  S,  ATMOSPHERIC 
ADMINISTRATION 
ENVIRONMENTAL  RESEARCH  LABS 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 

01CY  ATTN  R.  GRUBB 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 

P.O.  BOX  92957 

LOS  ANGELES,  CA  90009 


01CY 

ATTN 

I. 

GARFUNKEL 

01CY 

ATTN 

T. 

SALMI 

01CY 

ATTN 

V. 

JOSEPHSON 

01CY 

ATTN 

S. 

BOWER 

OlCY 

ATTN 

D. 

OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

01CY  ATTN  RADIO  SCIENCES 


LOS  ALAMOS  NATIONAL  LABORATORY 
MS  D438 

LOS  ALAMOS,  NM  87545 


OlCY 

ATTN 

S.P.  GARY 

OlCY 

ATTN 

J.  BOROVSKY 

SANDIA  LABORATORIES 

P.O.  BOX  5800 

ALBUQUERQUE,  NM 

87115 

OlCY 

ATTN 

V.  BROWN 

OlCY 

ATTN  ' 

A.  THORNBROUGH 

OlCY 

ATTN 

T.  WRIGHT 

OlCY 

ATTN 

D.  DABLGREN 

OlCY 

ATTN 

3141 

OlCY 

ATTN 

SPACE  PROJ  D1V 

SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  94550 

01CY  ATTN  B.  MURPHEY 

01CY  ATTN  T.  COOK 


AUSTIN  RESEARCH  ASSOCIATION,  INC. 
1901  RUTLAND  DRIVE 
AUSTIN,  TX  78758 

01CY  ATTN  L.  SLOAN 

01CY  ATTN  R.  THOMPSON 

BERKELEY  RESEARCH  ASSOCIATES.  INC. 


P.O.  BOX  983 
BERKELEY,  CA  94701 


OlCY 

ATTN 

J. 

WORKMAN 

OlCY 

ATTN 

C. 

PRETTIE 

OlCY 

ATTN 

s. 

BRECHT 

BOEING  COMPANY , 
P.O.  BOX  3707 

THE 

SEATTLE, 

WA  981 

24 

OlCY 

ATTN 

G. 

KEISTER 

OlCY 

ATTN 

D. 

MURRAY 

OlCY 

ATTN 

G. 

HALL 

OlCY 

ATTN 

J. 

KENNEY 

OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  DC  20545 

01CY  ATTN  DR.  YO  SONG 


CHARLES  STARK  DRAPER  LABORATORY 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  92139 

01CY  ATTN  D.B.  COX 

01CY  ATTN  J.P.  GILMORE 


COMSAT  LABORATORIES 
22300  COMSAT  DRIVE 
CLARKSBURG.  MD  20871 
01CY  ATTN  G.  HYDE 


36 


CORNELL  UNIVERSITY 

DEPT  OF  ELECTRICAL  ENGINEERING 

ITHACA,  NY  14850 

01CY.  ATTN  D.T.  FARLEY,  JR. 
ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 

01CY  ATTN  H.  LOGSTON 

OlCY  ATTN  SECURITY/ 

(PAUL  PHILLIPS) 

EOS  TECHNOLOGIES,  INC. 

606  WILSHIRE  BLVD. 

SANTA  MONICA,  CA  90401 

OlCY  ATTN  C.G.  GABBARD 

OlCY  ATTN  R.  LELEVIER 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

OlCY  ATTN  SECURITY  OFFICER 

OlCY  ATTN  T.N.  DAVIS 

OlCY  ATTN  NEAL  BROUN 

GTE  SYLVANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP- 
EASTERN  DIVISION 
77  A  STREET 
NEEDHAM,  MA  02194 

OlCY  ATTN  DICK  STEINHOF 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

OlCY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

OlCY  ATTN  DAN  MCCLELLAND 

OlCY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSIS 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA,  VA  22311 


OlCY 

ATTN 

J.M. 

AEIN 

OlCY 

ATTN 

ERNEST  BAUER 

OlCY 

ATTN 

HANS 

WOLFARD 

OlCY 

ATTN 

JOEL 

BENGSTON 

INTL  TELL  &  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
N'JTLEY,  NJ  07110 

OlCY  ATTN  TECHNICAL  LIBRARY 


JAYCOR 

P.O.  BOX  85154 
11011  TORREYANA  ROAD 
SAN  DIEGO,  CA  92138 

OlCY  ATTN  N.T.  GLADD 

OlCY  ATTN  J.L.  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

OlCY  ATTN  DOC  LIBRARIAN 

OlCY  ATTN  THOMAS  POTEMRA 

OlCY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORPORATION 
P.O.  BOX  7463 
COLORADO  SPRINGS,  CO  80933 
OlCY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED 
STUDIES 

816  STATE  STREET 
(P.O.  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 


OlCY 

ATTN 

DASIAC 

OlCY 

ATTN 

WARREN  S.  KNAPP 

OlCY 

ATTN 

WILLIAM  MCNAMARA 

OlCY 

ATTN 

B.  GAMBILL 

LINKABIT  CORPORATION 

10453  ROSELLE 

SAN  DIEGO,  CA  92121 

OlCY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 

OlCY  ATTN  DEPT  60-12 

OlCY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 


OlCY 

ATTN 

MARTIN  WALT/ 

DEPT  52-12 

OlCY 

ATTN 

W.L.  IMHOF/ 

DEPT.  52-12 

OlCY 

ATTN 

RICHARD  G.  JOHNSON/ 
DEPT.  52-12 

OlCY 

ATTN 

J.B.  CLADIS/ 

DEPT.  52-12 

MARTIN  MARIETTA  CORPORATION 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 

01 CY  ATTN  R.  HEFFNER 

MCDONNELL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  02647 


OlCY 

ATTN 

N.  HARRIS 

OlCY 

ATTN 

J.  MOULE 

OlCY 

ATTN 

GEORGE  MROZ 

OlCY 

ATTN 

V.  OLSON 

OlCY 

ATTN 

R.W.  HALPRIN 

OlCY 

ATTN 

TECHNICAL  LIBRARY 
SERVICES 

MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 


SANTA  BARBARA, 

CA  03101 

OlCY 

ATTN 

P.  FISCHER 

OlCY 

ATTN 

W.F.  CREVIER 

OlCY 

ATTN 

STEVEN  L.  GUTSCHE 

OlCY 

ATTN 

R.  BOGUSCH 

OlCY 

ATTN 

R.  HENDRICK 

OlCY 

ATTN 

RALPH  KILB 

OlCY 

ATTN 

DAVE  SOWLE 

OlCY 

ATTN 

F.  FAJEN 

OlCY 

ATTN 

M.  SCHEIBE 

OlCY 

ATTN 

CONRAD  L.  LONGMIRE 

OlCY 

ATTN 

B.  WHITE 

OlCY 

ATTN 

R.  STAGAT 

OlCY 

ATTN 

D.  KNEPP 

OlCY 

ATTN 

C.  RINO 

MISSION  RESEARCH  CORPORATION 
1720  RANDOLPH  ROAD,  S.E. 
ALBUQUERQUE,  NM  87106 

01CY  ATTN  R.  STELLINGVERF 

01CY  ATTN  M.  ALME 

01CY  ATTN  L.  WRIGHT 


PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD 
LOS  ANGELES,  CA  90025 

01CY  ATTN  E.C-  FIELD,  JR 
PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
UNIVERSITY  PARK,  PA  16802 
01  CY  ATTN  IONOSPHERIC 
RESEARCH  LAE 

PHOTOMETRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

OlCY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  04610 

OlCY  ATTN  A.  THOMSON 


R  &  D  ASSOCIATES 
P.O.  BOX  9695 


MARINA  DEL  REY, 

CA  90291 

OlCY 

ATTN 

FORREST  GILMORE 

OlCY 

ATTN 

W.B.  WRIGHT,  JR 

OlCY 

ATTN 

W.J.  KARZAS 

OlCY 

ATTN 

H.  ORY 

OlCY 

ATTN 

C.  MACDONALD 

OlCY 

ATTN 

BRIAN  LAMB 

OlCY 

ATTN 

MORGAN  GROVER 

RAYTHEON  CORPORATION 
528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

OlCY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42ND  STREET 
NEW  YORK,  NY  10036 

OlCY  ATTN  VINCE  TRAPANI 


MITRE  CORPORATION 
WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 

OlCY  ATTN  W.  HALL 

OlCY  ATTN  W.  FOSTER 


SCIENCE  APPLICATIONS 

INTERNATIONAL  CORPORATION 
10260  CAMPUS  POINT  DRIVE 
SAN  DIEGO,  CA  92121-1522 


OlCY 

ATTN 

L.M.  LINSON 

OlCY 

ATTN 

D. A.  HAMLIN 

OlCY 

ATTN 

E.  FRIEMAN 

OlCY 

ATTN 

E. A.  STRAKER 

OlCY 

ATTN 

C.A.  SMITH 
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SCIENCE  APPLICATIONS 

INTERNATIONAL  CORPORATION 
1710  GOODRIDGE  DRIVE 
MCLEAN,  VA  22102 

01CY  ATTN  T.  COCKAYNE 

01CY  ATTN  E.  HYMAN 


SRI  INTERNATIONAL 
333  RAVENSVOOD  AVENUE 


MENLO  PARK,  CA 

94025 

01CY 

ATTN 

J.  CASPER 

01CY 

ATTN 

DONALD  NEILSON 

01CY 

ATTN 

ALAN  BURNS 

01CY 

ATTN 

G.  SMITH 

01CY 

ATTN 

R.  TSUNODA 

01CY 

ATTN 

D.A.  JOHNSON 

01CY 

ATTN 

V.G.  CHESNUT 

01CY 

ATTN 

C.L.  RING 

01CY 

ATTN 

VALTER  JAYE 

01CY 

ATTN 

J.  VICKREY 

01CY 

ATTN 

R.L.  LEADABRAND 

01CY 

ATTN 

G.  CARPENTER 

01CY 

ATTN 

G.  PRICE 

01CY 

ATTN 

R.  LIVINGSTON 

01CY 

ATTN 

V.  GONZALES 

01CY 

ATTN 

D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  VIGGINS  AVENUE 
BEDFORD,  MA  01730 

01CY  ATTN  V.P.  BOQUIST 

TRW  DEFENSE  &  SPACE  SYS  GROUP 

ONE  SPACE  PARK 

REDONDO  BEACH,  CA  90278 

01CY  ATTN  R.K.  PLEBUCH 

01CY  ATTN  S.  ALTSCHULER 
01CY  ATTN  D.  DEE 
01CY  ATTN  D.  STOCKVELL/ 
SNTF/1575 

VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MA  01803 
01CY  ATTN  V.  REIDY 
01CY  ATTN  J.  CARPENTER 
01CY  ATTN  C.  HUMPHREY 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 

01CY  ATTN  N.  ZABUSKY 
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